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Abstract 

Magnetic monopoles have attracted the attention of physicists since the founding 
of the electromagnetic theory. Their search has been a constant endeavor which was 
intensified when Dirac established the relation between the existence of monopoles and 
charge quantization. However, these searches have been unsuccessful. We have recently 
proposed that monopolium, a monopole-antimonopole bound state, so strongly bound 
that it has a relatively small mass, could be easier to find and become an indirect but 
clear signature for the existence of magnetic monopoles. In here we extend our previous 
analysis for its production to two photon fusion at LHC energies. 
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1 Introduction 



The theoretical justification for the existence of classical magnetic poles is that they 
add symmetry to Maxwell's equations and explain charge quantization [H EJ [3] . Dirac 
formulated his theory of monopoles considering them point-like particles and quantum 
mechanical consistency conditions lead to the so called Dirac Quantization Condition 



where e is the electron charge, g the monopole magnetic charge and we use natural units 
h = c=l. 

Numerous experimental searches for monopoles have been carried out but all have met 
with failure [H El El El El El ED]- The last, carried out by the CDF collaboration at the 
Fermi National Laboratory [11], found no monopoles and established a lower mass limit 
of 360 GeV. 

This lack of experimental confirmation has led many physicists to abandon the hope 
in their existence. A way out of this impasse is the old idea of Dirac [H [12], namely, 
monopoles are not seen freely because they are confined by their strong magnetic forces 
forming a bound state called monopolium [131 [H]. This idea was the leitmotiv behind our 
recent research [15] , namely we proposed that monopolium might be easier to detect than 
free monopoles. We showed that certain parameterizations of the mass and the width, 
allowed for such a scenario. 

The Large Hadron Collider (LHC) will soon enter in operation and will probe the new 
energy frontier opening possibilities for new physics including the discovery of magnetic 
monopoles either directly, a possibility contemplated long time ago, [TBI [IT], or through the 
discovery of monopolium, as we have been advocating. This development motivates our 
present research which analyzes the production of monopolium at LHC by the mechanism 
of photon fusion. 

Recently, Dougall and Wick have presented a calculation of monopole-antimonopole 
production from photon fusion at proton colliders [HI [19] . The mass limit of CDF was 
obtained assuming Drell-Yan (DY) production which dominates over other processes at 
Tevatron energies. Dougall and Wick have shown that at LHC energies photon fusion is 
the dominant process [TEIES]- We proceed to calculate the production of monopolium by 
means of this process and compare our results with monopole-antimonopole production. 
We conclude that, if monopolium is a strongly bound state, its cross section is larger than 
that for creating a heavy monopole-antimonopole pair and consequently gives rise to a 
more clear experimental signal. 

2 Monopolium production 

A useful computational theory of monopoles does not currently exist to perform a direct 
production calculation. For this reason we will employ a minimal model of monopole inter- 
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Figure 1: Diagrammatic description of the reaction studied. 



action which assumes a monopole photon-coupling which is proportional to the monopole's 
induced electric field g/3 for a monopole moving with velocity j3 [TB I [T9 l [20] . This approxi- 
mation can be shown to be almost equivalent to the low energy effective theory of Ginzburg 
and Schiller [2TJ [22] . This theory was derived from the standard electroweak theory in the 
one loop approximation leading to an effective coupling proportional to g e ff ~ — g, where 
uj is a kinematical energy scale of the process which is below the monopole production 
threshold, thus rendering the theory perturbative. This is so because that in a photon 
fusion diagram the dynamical scale is \/ E 2 — Am 2 , thus 

uj VE 2 - Am 2 E(3 

— ~ ~ , (2) 

m 2m 2m 

and therefore if E ~ 2m, i.e., kinetic terms are small, both schemes coincide. Here uj 
describes the active energy scale, E the center of mass energy, m the monopole mass and 
(5 the velocity. 

The Dirac quantization condition does not specify the spin of the monopoles. We 
choose here monopoles of spin 1/2, following Dougall and Wick [181 [19J , coupled in mo- 
nopolium to spin in order to have a minimum energy radial structure. 

We then study the expected production process at LHC, namely 

p + p -> p(X)+p(X) + M, (3) 

shown in FigfTJ where p represents the proton, X an unknown final state and M the 
monopolium. This diagram summarizes the three possible processes: 

i) inelastic p + p -> X + X + (77) -> X + X + M 

ii) semi-elastic p + p^p + X + (77) — >• p + X + M 

iii) elastic p + p^p + p+ (77) — > p + p + M. 
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Figure 2: Diagrammatic description of the elementary subprocess of the monopolium production 
from photon fusion. 



In the inelastic scattering, both intermediate photons are radiated from partons (quarks 
or antiquarks) in the colliding protons. 

In the semi-elastic scattering one intermediate photon is radiated by a quark (or 
antiquark), as in the inelastic process, while the second photon is radiated from the 
other proton, coupling to the total proton charge and leaving a final state proton intact. 

In the elastic scattering both intermediate photons are radiated from the interacting 
protons leaving both protons intact in the final state. 

The full 77 calculation includes contributions from these three individual regimes. 

We next proceed to describe the elementary subprocess shown in Fig. El which deals 
only with photons and monopolium, and will come back to the full pp scattering treatment 
later on. The standard expression for the cross section of this elementary subprocess 
results in 

a(2l _ M ) = % M2r 'f» F " - (4) 
V ; E 2 (E 2 - M 2 ) + M 2 T 2 M y ' 

where we have assumed that monopolium decays with a width Tm and T(E), with E 
off mass shell, describes the production cross section. The width Tm arises from the 
softening of a delta function, 5(E 2 — M 2 ) and therefore is, in principle, independent of 
the production rate T(E) [23] . 

We enter now the computation of the T(E), which represents the width of the 27 
decay of monopolium. The calculation, following standard field theoretic techniques of 
the decay of a non-relativistic bound state [23J, EI] , leads to 

r(£) = ^hM0)| 2 . (5) 

We have used the conventional approximations, namely that the monopoles are almost on 
shell and that in the calculation of the elementary process we have neglected the binding 
energy, i.e. M = 2m. However, we express the final formula in terms of the monopolium 
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Figure 3: Diagramatic representation of the model describing the coupling of photons to mo 
nopolium. 



mass, M, because the latter will take care of the binding. Here a g corresponds to the 
photon-monopole coupling and ipu is the mo nopolium ground state wave function. 

Using the Coulomb wave functions of ref . [TB] expressed in the most convenient way to 
avoid details of the interaction, which will be parameterized by the binding energy, one 
has 



l<Mo)p 



- 2 3/2 m 3 

7r m 



(6) 



and the scheme of Dougall and Wick [TH [19] adapted to monopolium production, gives 
rise to 



r(E) 



(7) 



Here, a is the fine structure constant and /3 the monopolium velocity, 




0=M l-Tsr 



which is the velocity of the monopoles moving in the monopolium system. 

Note that due to the value of j3 the width vanishes at the monopolium mass, where 
the velocity is zero. Therefore a static monopolium is stable under this interaction. 

A caveat must be made. There is a duality of treatments in the above formulation, 
see Fig. El The static coupling is treated as a Coloumb like interaction of coupling g 
binding the monopoles into monopolium, although ultimately the details are eliminated 
in favor of the binding energy parameterized by the monopolium mass M. We find in this 
way a simple parametric description of the bound state. The dynamics of the production 
of the virtual monopoles, to be bound in monopolium, is described in accordance with 
the effective theory [T8l [19], and this coupling is f3g. This is similar to what is done 
in heavy quark physics [25j(see his figure 5), where the wave function is obtained by a 
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Figure 4: We show the photon fusion production cross section for monopolium (solid curve) 
for R = 1.5 and Tm = 0.1, and for monopole-antimonopole (dotted curve) as a function of the 
energy variable £ = E/M. 



parametric description using approximate strong dynamics while the coupling to photons 
is elementary. 

The production cross section can now be written as 

M a(27 ^ M) = ^ (e*-i)* + rv { ) 

where R = 2m/M, (1 < R < oo). This parameter ratio describes the binding energy of 
monopolium in units of M, since M = 2m + Ending. £ = E/M is the center of mass 
energy measured in units of M and T M = T M /M is the decay width of monopolium 
also measured in units of M. The right hand side is adimensional and therefore the 
above expression gives the cross section in units of 1/M 2 . The monopolium width, which 
vanishes with our dynamics, arises from higher order effects or other possible dynamics 
and we consider it as a parameter. Even in the case monopolium would be stable at rest 
the width would parameterize the beam width [23] . 

In Fig. H] we show the photon fusion cross section for monopolium for a value of 
R = 1.5 and r&r = 0.1, together with that for monopole-antimonopole obtained in ref. 
[HI [19], which we have simply rewritten using M as the energy unit, 

MV ^ - = log - (2 - " ' ■ (10) 



where ? = - 
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The qualitative features of both cross sections are radically different. The monopolium 
cross section is a spike slightly above the monopolium mass, where it vanishes, assuming a 
reasonable width Tm < M, while the monopole-antimonopole cross section is a soft curve 
extending over a large energy region. It is also clear from Eq. ([9]) that the height of the 
pick strongly depends on the value of the binding, i.e., the larger the binding energy the 
larger will be R and the higher will be the pick. 

It is important to note that the approach followed to describe the coupling (g — > g/3) 
provides us with a large negative power of energy which makes the monopolium cross 
section fall off very rapidly. In the Ginzburg- Schiller approach, four powers of the energy 
are substituted by four powers of the monopolium mass and therefore one obtains a larger 
effective width for the pick, and thus a larger integrated cross section for the same values 
of R and T M . 

3 Cross section estimates 

We calculate 77 fusion for monopolium production following the formalism of Drees et al. 
[26] benefitting from the the full documentation of the calculation in the work of Dougall 
and Wick [181 US]- We obtain therewith the pp cross section for monopolium and for 
monopole-antimonopole production within the same computational codes. 

The full pp calculation includes contributions of three types: inelastic, semi-elastic, 
and elastic scattering. We sum these individual contributions to find the total pp cross- 
section, o to t- 

In the inelastic scattering, p + p — > X + X + (77) — > X + X + M, to approximate the 
quark distribution within the proton we use the Cteq6-1L parton distribution functions 
[27] and choose Q 2 = s/4 throughout. 

We employ an equivalent-photon approximation for the photon spectrum of the in- 
termediate quarks [2"gj |2"§]. 

In semi-elastic scattering, p + p^p + X+ (77) — > p + X + M, the photon spectrum 
associated with the interacting proton must be altered from the equivalent-photon ap- 
proximation for quarks to account for the proton structure. To accommodate the proton 
structure we use the modified equivalent-photon approximation of [30J. 

For elastic scattering, p + p — > p + p+ (77) — > p + p + M, both protons remain intact 
in the final state. 

In Fig{5] we show the comparison between monopole-antimonopole and monopolium 
production cross sections using M as energy unit as a function of R = 2m/ M, i.e., binding 
energy \E binding \/M = R — 1. Two regimes are found: 

i) the low binding regime 1 < R < 2 where the mm production is dominant; 

ii) the large binding regime R > 2 where monopolium production is dominant and can 
be very large. 

In Figj6] we plot the logarithm of the total cross section (in fb) for pp production of 
monopolium as a function of monopole mass, for two masses of monopolium (M=100 
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Figure 5: The left figure shows logarithmic plots of the pp production cross sections (in units 
of M 2 ) via photon fusion for monopolium (Tm = 0.1M) and for monopole-antimonopole as a 
function R = 2m /M for low binding. The figure on the right shows logarithmic plots of the 
pp production cross section (in units of M 2 ) via photon fusion for monopolium (Tm = 0.1M) 
and for monopole-antimonopole as a function of R for large binding. The solid line represents 
the total cross section; the various contributions are represented by: inelastic (dashed line); 
semielastic (dotted line) and elastic (dot-dashed line). 

GeV and M= 1000 GeV), as well as the total cross section for monopole-antimonopole 
production as a function of monopole mass. The two mass plot shows the dependence 
with monopolium mass of the cross section. There is always a threshold at m = M/2, 
i.e. for zero binding energy, where the cross section vanishes, and then a rapid rise with 
monopole mass. Thereafter the cross section grows in a softer manner, although it should 
be realized that we deal with a log plot and therefore the growth is not negligible. Finally, 
the smaller the monopolium mass, the larger the cross section is. 

The most spectacular signature of the curves is that for fixed monopolium mass 
the cross section increases with monopole mass, instead of decreasing, as happens in 
monopole-antimonopole production, and the magnitude becomes enormous for very large 
binding energy. Thus a strongly bound monopolium state would be an ideal system to 
disentangle monopole dynamics. 

4 Conclusions 

We have carried out an investigation looking for hints of the so far not seen monopoles. 
Our motivation has been, that of our previous work [15] . namely that monopolium, if 
strongly bound, is easier to produce than monopole-antimonopole pairs. 

We have performed a calculation in which monopolium is produced via the conven- 
tional monopole dynamics [HI [19] used to study monopole-antimonopole production by 
photon fusion at LHC energies. Our calculation is parameterized in terms of three quan- 
tities: m, the monopole mass, M the monopolium mass (or binding energy) and the width 



7 



15 




500 1000 1500 2000 2500 3000 

m (GeV) 

Figure 6: We represent the logarithmic plot of the total pp cross section from photon fusion 
in femtobarns as a function of monopole mass, for two monopolium masses M = 100 GeV and 
M = 1000 GeV and compare it with the corresponding monopole-antimonopole (m m) cross 
section. 

r^f which we simply take to be narrow compared with the monopolium mass. 

Our analysis distinguishes two clearly distinct regions associated with the binding 
energy in monopolium. If the binding energy is small compared to the monopole mass, 
the monopole-antimonopole process will be dominant. On the contrary if it is comparable 
or larger, the monopolium process is not only dominant but can be extremely large. 

The most favorable scenario, which agrees with that discussed in our previous work, 
is a two energy scale scenario, whose 

i) low energy scale is governed by the monopolium mass, M, reachable by LHC, 

ii) and whose high energy scale is governed by the monopole mass, m, which arises 
through the structure of monopolium, and which could be larger than the energy 
reachable by LHC. 

Under these circumstances the cross section as a function of the monopole mass be- 
comes sizeable. 

An example, if monopole and monopolium have a mass of 1 TeV, at an integrated 
luminosity of 100 fb^ 1 at LHC, Dougall and Wick predict about 10 6 monopoles, while 
our calculation would produce 10 8 monopolia. 

Since at present we cannot calculate the monopolium parameters, M and Tm, the 
experimental endeavor is not easy. However, if we extend the theory to incorporate the 
weak interaction [2TJ [22] there are some features which might simplify the task, 
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i) the resonance peak of the monopolium can be found in three exit channels 27, 7 Zq 
and 2 Zqs. 

ii) monopolium can be produced in an excited state before it annihilates, thus the 
annihilation process will be accompanied by a Rydberg radiation spectrum; 

The calculated values for the cross sections, corresponding to reasonable monopolium 
mass scenarios, render our calculation sound and this line of research worth pursuing. 
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